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Abstract. Under conditions of hmited fixed-nitrogen, some filamentous cyanobac- 
teria develop a regular pattern of heterocyst cells that fix nitrogen for the remaining 
vegetative cells. We examine three different heterocyst placement strategies by quanti- 
tatively modelling filament growth while varying both external fixed-nitrogen and leak- 
age from the filament. We find that there is an optimum heterocyst frequency which 
maximizes the growth rate of the filament; the optimum frequency decreases as the ex- 
ternal fixed-nitrogen concentration increases but increases as the leakage increases. In 
the presence of leakage, filaments implementing a local heterocyst placement strategy 
grow significantly faster than filaments implementing random heterocyst placement 
strategies. With no extracellular fixed-nitrogen, consistent with recent experimental 
studies of Anabaena sp. PCC 7120, the modelled heterocyst spacing distribution us- 
ing our local heterocyst placement strategy is qualitatively similar to experimentally 
observed patterns. As external fixed-nitrogen is increased, the spacing distribution for 
our local placement strategy retains the same shape while the average spacing between 
heterocysts continuously increases. 
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1. Introduction 

Cyanobacteria are prokaryotes able to grow photoautotrophically using oxygenic 
photosynthesis. They preferentially use ammonium or nitrate as sources of fixed- 
nitrogen [l]. However, in conditions of low exogenous fixed-nitrogen, cyanobacteria can 
fix atmospheric nitrogen. Heterocystous filamentous cyanobacteria form unbranched 
clonal filaments of cells, and fix nitrogen within dedicated and terminally differentiated 
heterocyst cells [21 [3] . Heterocysts are separated by clusters of photosynthetic vegetative 
cells that cannot themselves produce fixed-nitrogen (fN). Significant progress has 
been made in characterizing the genetic network underlying the observed heterocyst 
pattern [21 H] that is developed in response to the absence of external fixed-nitrogen. 
Nevertheless, little attention has been paid to the functional role of the heterocyst 
pattern itself in these model developmental organisms. 

As we shall show, a filament with too few (or inactive) heterocysts will starve 
of fixed-nitrogen and grow slowly without external fixed-nitrogen [5J while too many 
heterocysts, which do not grow and divide, will also inhibit growth [6]. Balancing these 
effects in the filament qualitatively explains the heterocyst frequency of approximately 
10% observed with no external fixed-nitrogen, but does not explain the distinctive 
pattern of heterocyst spacings that is seen in the model cyanobacterium Anabaena sp. 
PCC 7120 [HI [71 El E] . Indeed, the observation that mutant strains exhibiting a distinct 
pattern of multiple- contiguous heterocysts (Mch) show reduced growth [10] indicates 
that heterocyst placement is important. In this paper, we use quantitative modelling 
to explore the hypotheses that simple heterocyst placement strategies can affect 
filament growth, and that observed heterocyst patterns reflect placement strategies that 
maximize growth. 

Previous models of cyanobacterial spacings have directly compared model and 
experimental distributions. Meeks and Elhai [H] compared the expected heterocyst 
spacing for randomly spaced heterocysts to the experimentally measured distribution 
and found clear disagreement. Wo Ik and Quine [12] examined a diffusible inhibitor 
mechanism in which a radius of inhibition around existing heterocysts depends on 
diffusion and degradation of the inhibitor, and obtained qualitative agreement with 
early observed spacing distributions for Anabaena sp. PCC 7120 (hereafter simply PCC 
7120). These models identify lateral inhibition as a plausible mechanism behind the 
observed heterocyst patterns but could not examine why (or if) the observed pattern is 
optimal, or how the pattern may change under different experimental conditions. 

We hypothesize that leakage of fN from the cyanobacterial filament may distinguish 
between heterocyst patterns, in terms of growth rates. Evidence for leakage was 
first reported by Fogg et al [131 HI] 5 who found fN products outside the filament. 
Supporting this, Paerl [15] observed bacteria clustering around cyanobacterial filaments, 
particularly at junctions between heterocysts and vegetative cells, and indicated a 
possible symbiosis based on enhanced nitrogen fixation and leaked fixation products. 
Thiel [16] found protein proteolysis byproducts in the extracellular medium after fixed- 
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nitrogen starvation. Significant leakage is also consistent with the observation of reduced 
diazotrophic growth in strains with impaired amino acid uptake transporters pTl ITS]. 
Nevertheless, we are not aware of any quantitative measurements of the rate of leakage. 
Qualitatively, a regular pattern of heterocysts would minimize the distance traveled by 
fN products so as to minimize leakage from vegetative cells, leaving more fN available 
for growth. 

In this paper, we examine heterocyst frequencies that maximize filament growth 
within the context of a quantitative transport model [19] that incorporates fixed-nitrogen 
transport, vegetative cell growth, and fixed-nitrogen production at heterocysts. We 
explore the impact on growth of different heterocyst placement strategies, including 
random placement, and find that they are almost indistinguishable in the absence of 
leakage of fN from the filament, but clearly distinct with leakage. We find that the 
heterocyst spacing patterns corresponding to maximal filament growth are qualitatively 
similar to those seen experimentally. 

The genera Anabaena, which includes heterocystous cyanobacteria, is widely 
distributed geographically in freshwater lakes [201 121] where fixed-nitrogen is one of 
the major substances limiting growth [22l [23l El]. Nitrification turns ammonia into 
nitrite and then nitrate [21], [25]. Nitrate levels in lakes range from 0-10 mg/1 (0- 
10x10^^ m~^) in unpolluted freshwater, but vary both seasonally and spatially [2T| [26]. 
Heterocystous cyanobacteria are also found in the oceans [271 [28] where nitrate levels 
have been recorded to vary from approximately 0-50 /xmol/kg [29] (0-3x10^^ m~^). 
It has long been known that sources of fixed-nitrogen initially present in the medium 
can increase the mean spacing between heterocysts [131 SO] and that heterocysts will 
differentiate at non-zero levels of external fixed- nitrogen (efN) [13]. Field studies of 
heterocystous cyanobacteria also show both a range of heterocyst counts and a range 
of efN concentrations, with significant positive correlations between the two (see e.g. 
[31]). In addition, steady-state chemostat experiments show that Anabaena flos- aquae 
can adjust nitrogen fixation to achieve constant growth, within approximately 10%, as 
efN concentrations are varied [32]. While growth independence from efN seems desirable 
in the face of environmental variability, it raises the question of how and how well it is 
achieved in terms of heterocyst fraction and pattern. 

The heterocyst pattern changes with time as it evolves towards a steady-state 
distribution after efN deprivation [HI [H [9]. Earlier models [HI [12] have focused on 
the early heterocyst pattern observed 24 h after efN deprivation. While the early 
pattern and the later patterns are qualitatively similar, with broad distributions of 
heterocyst spacings ranging from zero to more than twenty cells between heterocysts 
with a peak at about a ten cell spacing, we focus on the steady-state distribution in this 
paper. We explore the hypothesis that cyanobacterial filaments use a unified heterocyst 
placement strategy even at non-zero levels of efN. With all of our strategies, we find 
that maximal growth is observed with heterocyst frequencies that decrease continuously 
with increasing levels of efN. 
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2. Model 

2.1. Fixed-nitrogen Transport 

Our model of fixed-nitrogen (fN) transport and incorporation is adapted from Brown 
and Rutenberg [19], with the addition of import from the extracellular medium. That 
study found that periplasmic transport was not required to explain nanoSlMS studies 
of fixed-nitrogen profiles [33], so we do not include periplasmic transport in our model. 
Our model tracks the total amount of freely-diffusing fN, N{i, t), in each cell i vs. time 
t: 

j^N{i, t) = ^tot{i, t) + Djp,f,,l{i, t) - DlN{i, t) + G„ (1) 

where ^totii) = - 1) + ^hi^ + 1) - ^hij) - is the net diffusive flux into cell 

i, equal to the sum of the incoming fluxes minus the outgoing fluxes. Each cell has 
two outgoing fluxes, $l and to its left and right neighbours, respectively, and two 
incoming fluxes from its neighbours. Following Pick's law, each flux is the product of 
the local density N{i,t)/l{i,t) and a transport coefficient Dc- We use Dc=l.54: fim 
s"^ between two vegetative cells and Dc=0.19 fim s~^ between a vegetative cell and a 
heterocyst [I9] — we keep these values fixed in this paper. 




Figure 1. Schematic of fixed-nitrogen transport as represented by Eqn. [T] N{i) is 
the amount of freely-diffusing fixed-nitrogen for ceh i. Dc governs transport between 
cehs, Dj governs import from the external medium into the cell, and Dl controls loss 
from the cell to the external medium. 

In addition to fluxes along the filament, Dj is the coefficient for import from the 
external medium with efN concentration Pe/Af to the cell of length /j and -Dl is the 
coefficient for losses from the cell to outside the filament. These transport processes are 
shown schematically in Fig. [H 

To estimate Dj, we use Fogg's observations of heterocyst frequency in a medium 
to which a fixed amount of efN in the form of ammonia was added [13]. Heterocyst 
frequency decreased in time due to vegetative growth until the ammonia concentration 
dropped to 4x 10~^M (see p. 245 of [13]). We assume that a slightly higher concentration 
of efN, 5 X 10"^M (pe/Af=3xl0^^m~^), is just sufficient for maximal growth. From 
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Eqn. [T] the amount of fixed-nitrogen imported from the external medium in a time 
T is DjpefN^T. For a cell to double in length over a period of time T=20h it needs 
approximately 1.4x10^° N atoms to be imported [I9]. Using /=/mm=2.25/xm, this yields 
Di = 2.9 X 10" V/(/im-s). 

We also include a loss term with coefficient -Dl in Eq. [T], following [131 El HSl [lEl 
[TTllIH]. Loss is expected to be proportional to the cytoplasmic density, N/{Al), where 
I is the length of the cell and A is the cross-sectional area, and also to the number of 
transporters, which will themselves be proportional to the cell length. This results in a 
loss term that is independent of length /, with a coefficient with units of s^^. There 
are no direct measurements of leakage rates of fN from cyanobacterial filaments that 
we are aware of. We take the loss rate as small compared to the import rate, and so 
use either Dl=0.01Dj/A (referred to as 1% loss) and Dl=0.1Dj/A (referred to as 10% 
loss). We use a fixed A = vr/xm^ corresponding to a radius of l/im. 

In addition to transport terms, there is also a source/sink term Gi in Eqn. [1] 
that describes fN production and consumption in the heterocysts and vegetative cells, 
respectively. This G term is discussed in the next section. 

2.2. Cell Growth and Division 

Following [19], we take PCC 7120 cells to have a minimum size of Imin = 2.25/im and 
a maximum size of Imax = '^Imin [H E]. When a cell reaches Imax it is divided into two 
cells of equal length, each of which is randomly assigned a new growth rate and half of 
the fixed-nitrogen in the parent cell. We initialized lengths randomly from an analytical 
steady-state distribution of cell lengths ranging between /^m and Imax [M] and used an 
average doubling time To = 20h [T8l [35] . We define a minimum doubling time T^m = 
Td — A and a maximum doubling time T^ax = Tr, -(- A. A doubling time T is randomly 
and uniformly selected from this range and converted to a growth rate R=lmin/T . We 
take A = 4.5h [191 [36]. 

In Eqn. [H for heterocysts Gi = Ghet- The heterocyst fN production rate, 
Ghet = 3.15 X 10^s~\ is chosen to supply the growth of approximately 20 vegetative 
cells. 

In Eqn. [H for vegetative cells Gi = G^eg is a sink term determining the incorporation 
rate of cytoplasmic fN removed to support cellular growth. G^eg depends on the actual 
growth rate R of the cell (in /xm/s), which in turn depends upon the locally available 
cytoplasmic fN: 



where g is the amount of fixed-nitrogen a cell needs to grow per fim; g = 1.4x 10^°//mjn — 
6.2 X 10^/xm~^ [19]. We assume that cells increase their length at their maximal growth 
rate as long as there is available cytoplasmic fixed-nitrogen. Otherwise, they can 
only grow using the fixed-nitrogen fiux into the cell from neighbouring cells: 



G, 



veg 



gR{R^^\N{i,t)) 



(2) 




(3) 
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Note that cells with = may still grow, but will be limited by the incoming fluxes 
of fixed-nitrogen from both adjoining cells and the external medium, Gyeg = ^in = 
^R{i - 1) + + 1) + DjpefNl{i). 

2.3. Heterocyst Placement 

While we know that experimentally observed heterocyst patterns are not due to random 
placement [H], it is useful to evaluate the steady-state growth rate achievable with 
random heterocyst placement. We will use random placement as a point of reference 
for other heterocyst placement strategies. We consider three simple strategies: random 
placement, random placement with no contiguous heterocysts, and local placement. For 
all of them, once a heterocyst is placed it immediately stops growing {R"^^ = 0) and 
fixes nitrogen (Gj = Ghet)- Commitment to differentiation does not occur until after 
approximately 8 h of efN deprivation and can take as long as 14 h [H [371 [38]. Heterocysts 
mature and begin to fix nitrogen approximately 18-24 h after efN deprivation [3l M, [37] 
(though see [36j)- For random strategies, heterocyst placement corresponds to when 
heterocysts begin to produce fN. For local placement, we have an explicit delay (see 
below) . 

Our first, reference, strategy is random placement ("random"). The heterocyst 
fraction / is fixed, and during filament growth a random vegetative cell is replaced with 
a heterocyst whenever possible but without exceeding /. 

Our second strategy ("no-Mch") refiects the observation that contiguous 
heterocysts are not observed during normal development [2[ [3]. It consists of our 
random strategy, but with the additional restriction that vegetative cells adjacent to 
existing heterocysts are never selected for development. 

Our third strategy is local heterocyst placement ("local"). Any vegetative cell that 
has been continuously starving for a defined interval r is changed into a heterocyst. 
Starvation is defined as N{i) = 0, so that these cells have reduced growth R < for 
a significant period of time. We vary r within the range 1 — 20 h. Starvation occurs 
due to distance from heterocysts [19], but can also refiect local clusters of fast-growing 
cells. 

2.4. Some numerical details 

Periodic boundary conditions were used to minimize end effects. Filaments were 
initiated with no heterocysts and the different strategies were followed until a steady- 
state was reached. For the local strategy, a "no-Mch" rule was followed for the first 24 
hours of differentiation to reduce initial transients. 

The growth rate constant n was calculated every six simulated hours using the 
total length of the filament L{t), where L{t) = L{t — Qh)e^'^^. Heterocyst frequency 
was sampled every six simulated hours by dividing the number of heterocysts by 
the total number of cells in the filament. Heterocyst spacings were also recorded 
every six simulated hours. Measurements were averaged daily, and the results of ten 
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independently seeded runs were used to determine an overall daily average and standard 
deviation. All data shown is for the fifth day, which exhibits steady-state for the 
parameters explored (in comparison with data from the fourth day). 

We begin each simulation with 100 cells. Most of the phenomena we investigate 
occurs in filaments with a growth rate constant greater than /i =0.02 h~^, which 
would allow the filament to grow to more than 1100 cells after five days. We use a 
computational timestep At =0.01s with a simpler Euler discretization of Eqn. [U smaller 
timesteps yield indistinguishable results. 

For random heterocyst placement strategies, optimal heterocyst frequencies for 
growth without leakage in Fig. [2](c), and the corresponding growth rates in Fig. M^d), 
were found by locally using the standard Marquardt-Levenberg fit algorithm with a 
quadratic function to the left of the optimum frequency and a linear function to the 
right, with the two functions meeting at the optimum frequency. Optimum heterocyst 
frequencies with leakage in Figs. Hl^b) and (d), and the corresponding growth rates in 
Figs. 111(a) and (c), were found by a least squares quadratic fit near the maximal growth 
rate. 



3. Results 

3.1. No Leakage 

We first examined systems with randomly placed heterocysts in filaments with zero 
leakage. Qualitatively there are two growth regimes: starving or excess fixed-nitrogen. 
At small heterocyst fractions and small levels of efN, starved growth will be determined 
by the amount of fN produced by heterocysts as well as fN imported from outside the 
filament. The fN from heterocysts is proportional to their fraction /, while the fN 
imported from outside the filament is proportional to PejN- 

Ghet r , Di , , 

t'-starve — J H PefNi yi) 

9 ' ''avg 9 

where lavg is the length of an average cell and g is the amount of fN per unit length 
needed for growth. At large heterocyst fractions or with large levels of efN, sufficient 
fN is present for growth but only the vegetative cells, with fraction 1 — /, will grow: 

f-^excess j \^ J ji \p) 



''mm 



where the /n(2) factor is needed to convert the growth rate i? of a single cell to a 
growth rate exponent /i. Both of these limiting behaviors ^starve and ^excess are shown 
as straight black lines in Figs. |2]^a) and (b). 

The modelling results for filaments with different levels of efN and with zero leakage 
[Dl = 0), shown in Fig.|2](a), agree well with the limiting behaviours for small or large /. 
There is a clear optimal heterocyst frequency, /*, at which growth is maximal. Increasing 
the amount of efN shifts ^starve up, while ^excess is unaffected. As efN increases, the 
optimal frequency decreases and the corresponding maximal growth rate n* increases. 
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Figure 2. (a) growth rate exponent /i vs. heterocyst frequency / for systems 
with zero leak;age (D^ — 0) and random heterocyst placement. Red circles indicate 
PefN — 0, blue squares indicate pefN = 0.5 x lO^^m^^, and green diamonds indicate 
PefN = 1 X lO^^m"^. Black lines indicate asymptotic behavior from Eqns.|3]and[S] (b) 
highlights the data near the maximal /z in (a). The red dashed curve indicates the best 
fit near the indicated optimal heterocyst fraction /*. (c) optimum heterocyst frequency 
/* vs. PefN using random heterocyst placement with zero leakage. Numerical data is 
shown as red circles while the frequency described by the intersection of Eqns. |4]and 
[5] is shown as a solid black line, (d) growth rate /i vs. pefN- Red circles with a solid 
line indicate p,* from random heterocyst placement with zero leakage. Orange squares 
with a dashed orange line indicate /i for a filament with no heterocysts, while the 
black diamonds and black dashed line indicate p for a filament maintaining a constant 
heterocyst frequency that is optimal at PefN = 0. 

The modelled growth rates are somewhat below the limiting regimes near /*, as 
highlighted in Fig. Et^b). This growth deficit appears to be due to finite diffusivity of 
fN within the filament, which limits the reach of excess fN from regions with excess 
heterocysts in the face of continued expansion of vegetative regions due to ongoing 
growth. Indeed, the growth deficit disappears when Dc — )■ oo (data not shown). The 
growth deficit is also smaller when heterocysts are placed close to starving cells (see 
below). We do not yet have an analytical treatment of this growth deficit, though it is 
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intriguing. As a result of the growth deficit, the optimal heterocyst frequency from our 
quantitative model, /*, is larger than given by the intersection of fistm-ve and fiexcess, as 
shown in Fig. [2](c). 

The maximal growth rates, fi*, corresponding to the optimal heterocyst frequencies 
/*, are plotted with solid red hues vs. Pe/7V in Fig. M^d). The dashed orange lines 
show the growth rate with no heterocysts (/ = 0), which exhibits sharply reduced 
growth at smaller Pe/N- The dashed black line shows the growth rate exhibited by the 
fixed heterocyst fraction /o that is optimal for Pe/TV = 0, which exhibits a constant 
but reduced growth when Pe/Af > 0. In general, filaments that maintain an optimal 
heterocyst fraction by maintaining / = /* as pefN varies will outgrow filaments with 
any given fixed heterocyst fraction. 

3.2. Leakage 




f 



Figure 3. Growth rate constant vs. heterocyst frequency / for different heterocyst 
placement strategies and leakage levels with zero external fixed-nitrogen. Zero leakage 
is shown in red, 1% leakage is shown in blue, and f 0% leakage is shown in green. Solid 
lines are random heterocyst placement, dashed lines are no-Mch heterocyst placement. 
Red circles, blue squares, and green diamonds show data for local heterocyst placement 
for zero, 1%, and 10% leakage respectively. The local heterocyst placement data points 
show variation of the period of starvation until commitment r from 1-20 h, with shorter 
T exhibiting larger fj,. 

As shown in Fig. [3l with no leakage the random (solid red line) and no-Mch (dashed 
red line) heterocyst placement strategies have almost indistinguishable growth rates. 
Local heterocyst placement (red circles) leads to slightly faster growth with short wait 
times r = 1 h, but shows significantly slower growth at longer wait times comparable 
to heterocyst maturation times. Note that / is not independently controlled with local 
placement, so that only a narrow range of / is seen as r is varied. Faster local growth is 
seen with smaller r as starving cells are provided with fN earlier, and experience shorter 
periods of restricted growth. 
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With non-zero leakage, the three heterocyst placement strategies produce 
noticeably different growth curves. At 1% leakage (blue lines and squares), the growth 
rates with no-Mch are slightly above the random placement strategy. More striking is 
the dramatic improvement with local placement, which has significantly better growth 
at any given heterocyst frequency, but also better maximal growth n* for r e [1, 13]h. 
The same trends continue when leakage is increased to 10% (green lines and green 
diamonds): the growth rates of filaments with local heterocyst placement are higher 
than the maximal growth rates of other strategies for r G [1, 17]h, and we also see that 
the corresponding heterocyst frequencies /* are much lower with local strategies than 
with either random placement or no-Mch placement strategies. 

Note that /* and fi* correspond to the heterocyst frequency and growth rate 
exhibited by the local placement strategy for a given delay time r, or to the growth- 
optimized heterocyst frequency and corresponding growth rate exhibited by the random 
or no-Mch placement strategies. 

In general, leakage of fN from the cyanobacterial filament will more strongly inhibit 
growth if heterocysts are not placed close to starving vegetative cells. In any case, 
leakage decreases the growth rate at a given heterocyst frequency, and so leads to a 
larger heterocyst frequency /* and correspondingly decreased growth rate fi* with both 
the local strategy and the growth-optimized random strategies. 

3.3. Varying External Fixed-nitrogen Concentration 

When the external fixed-nitrogen concentration, Pe/N, is increased from zero the 
differences between the heterocyst placement strategies diminish. Fig. |l](a) shows 
growth, fi*, vs. pefN for all three placement strategies with 1% leakage. Also shown 
with the orange-dashed line is the growth expected with no heterocysts (/ = 0). We see 
that local heterocyst placement has a significantly higher growth rate than optimized 
random or no-Mch strategies. The differences are largest at Pe/Af = 0, decrease as efN 
increases, and vanish when all growth is supported by efN alone. Similar qualitative 
behavior is seen with 10% leakage, in Fig. Hl^c), though the growth rates are somewhat 
lower. 

Figs. 111(b) and (d) show the heterocyst frequency, /*, for all three heterocyst 
placement strategies for a range of efN, for 1% and 10% leakage respectively. For most 
of the lower fixed-nitrogen levels, local heterocyst placement has a significantly lower 
heterocyst frequency than optimized random or no-Mch heterocyst placement. When 
heterocysts are placed near starving cells fewer heterocysts are necessary to satisfy the 
fixed-nitrogen requirements for the growing vegetative cells. As /* approaches zero, the 
trend reverses. The local heterocyst placement has a higher /* than either random or no- 
Mch placement because local can continue to beneficially place heterocysts near starving 
cells as / increases and fewer cells starve, while random strategies cannot. The reversal of 
this trend is particularly apparent with larger leakage levels (Fig.|4](d)), where placement 
of heterocysts close to fast-growing vegetative cells is particularly important, and leads 
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Figure 4. (a) growth rate constant fi vs. external fixed-nitrogen concentration pefN 
for filaments with 1% leakage. Red circles with a solid line are random heterocyst 
placement at the optimum frequency, blue squares with a solid line are no-Mch 
heterocyst placement at the optimum frequency, green diamonds with a solid line are 
local heterocyst placement, and the orange dotted line is filaments with no heterocysts. 
(b) heterocyst frequency at which maximum growth occurs vs. external fixed-nitrogen 
concentration PefN- Red circles with solid line are random heterocyst placement, blue 
squares with solid line are no-Mch heterocyst placement and green diamonds with solid 
line are local heterocyst placement, (c) and (d) are similar runs labeled in the same 
fashion as (a) and (b), respectively, with higher leakage of 10%. 



to a noticeable non-linearity of the heterocyst frequency /* with local placement at 
larger p^jN- 

3-4- Heterocyst Spacing 

Figs. [5]^a)-(c) show heterocyst spacing distributions for the different heterocyst 
placement strategies with 1% leakage, all with approximately / = 0.1 for ease of 
comparison. Both random placement and no-Mch favour small heterocyst spacings. 
Random placement, in Fig. |5t^a), peaks at adjacent heterocysts (corresponding to the 
Mch phenotype) and drops off for larger separations, while no-Mch placement, in 
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Fig.[S](b), peaks at a spacing of 4 intercalating cells between heterocysts. Significant bias 
towards even spacings, due to ongoing filament growth, is also seen. The distribution 
for local placement is quite different, with a symmetric peak at approximately 12 
intercalating cells and very few heterocysts separated by less than 6 cells. Fig. [5]^d) 
is an experimental spacing distribution after 96h of fixed-nitrogen deprivation for WT 
Anabaena PCC 7120 |8]. The heterocyst spacing distribution with the local placement 
strategy is qualitatively similar to the experimental distribution, even though the 
strategy selection was done with respect to growth alone — without consideration of 
the spacing distribution. 



(a) 




(b) 



Heterocyst Separation 




Heterocyst Separation 




Figure 5. (a)-(c) are heterocyst spacing distributions from stochastic simulations 
with 1% leakage and zero external fixed-nitrogen with error bars representing the 
standard deviation of the mean, (a) is random heterocyst placement with 10% 
heterocysts, (b) is no-Mch heterocyst placement with 10% heterocysts, and (c) is local 
heterocyst placement with period of starvation until commitment of r = 8h. (d) is the 
experimental steady-state WT distribution from [8] . 

With non-zero Pe/Af we consider only the heterocyst spacing distributions for the 
local heterocyst placement strategy, as illustrated in Figs. |6t^a)-(c) with 1% leakage. 
The peak separation increases with PefN, consistent with the decreasing /* we saw in 
Fig. m (b). The distribution also becomes significantly wider. By scaling the heterocyst 
separation by the average separation for each distribution we see in Fig. [6] (d) that 
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the distributions at different pefN approximately collapse to a single scaled distribution 
independent of PefN- The inset in Fig. |6t^d) is the average heterocyst separation vs. 

PefN- 
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Figure 6. Heterocyst spaeing distributions for filaments with local heterocyst 
placement, 1% leakage for various PefN and a period of starvation until committment 
of T = 8h. (a) PefN = 0.5 X lO^^m-^, (b) pefw = 1 x lO^^m-^^ and (c) 
PefN = 1.5 X lO^^m"'^. (d) shows the probability density of the distributions with 
PefN = (red), 0.5 (blue), 1 (green), and 1.5 (orange) xlO^^m"'^ vs. scaled heterocyst 
separation. Inset in (d) is the average heterocyst separation vs. efN concentration. 



4. Discussion 



We have explored three heterocyst placement strategies in a cyanobacterial filament 
with a quantitative model for fixed-nitrogen (fN) transport and dynamics that includes 
growth of vegetative cells and production of fN by heterocysts. The three strategies 
were random placement, random placement with no-Mch, and local placement. For 
random strategies, we found an optimal heterocyst frequency /* at which the growth 
rate of the filament is maximized, fi*. With fewer heterocysts than /* filaments are 
fN limited. With more heterocysts than /* the filament is hampered by an excess of 
non-growing heterocysts. As the external fN, pefN, is increased from zero, the optimal 
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heterocyst frequency decreases continuously until it reaches zero at a concentration p*j7v 
where all fN needs are met by imported extracellular fN. The local placement strategy 
led to a similar decreasing heterocyst fraction /* with PefN, though without explicit 
optimization of the heterocyst frequency. 

Without leakage of cytoplasmic fN from the filament, via Dl, the difference between 
these strategies is small (see Fig. [3], red lines and circles). With leakage there are 
significant growth differences between the strategies under limiting efN conditions 
(with PefN < pIjn)- Filaments using local placement of heterocysts grow faster than 
filaments with no-Mch, which in turn grow faster than filaments using random heterocyst 
placement. The differences between the strategies are largest when pefN = 0, and 
decrease as Pe/Af increases — until it vanishes above p*j7v where the heterocyst frequency 
/* = for all three strategies. We suggest that leakage of fN from growing filaments may 
be important for understanding the adaptive nature of heterocyst placement strategies. 

We beheve that the observed growth differences are relevant. Selective sweeps occur 
when a beneficial mutation 'sweeps' a population and becomes fixed [39] . Selective 
sweeps have been observed in cultures of E. coli for mutations with a fitness parameter 
as small as m = 0.006 [39], which in our system corresponds to a growth rate constant 
difference of 0.0003h~^. Fig. HJ^a) shows that the growth rate difference between random 
heterocyst placement strategies and local heterocyst placement is about 0.0025h~^ 
for 1% leakage for zero external fixed-nitrogen, and is similar for 10% leakage — 
substantially larger than necessary for a selective sweep in E. coli. 

Significantly, we did not pick our best strategy to match the heterocyst spacing 
pattern nor did we tweak the heterocyst frequency /. Rather we implemented a 
simple local strategy which itself chose /*. This best, local strategy resulted in 
faster growth than optimized random strategies, and in heterocyst spacings remarkably 
similar to those seen experimentally. The similarity of our model heterocyst spacing 
distributions with observed patterns suggests that local fN starvation may drive 
heterocyst development. Indeed, we believe that the extensive genetic network of 
Anabaena heterocyst development [3] effectively implements something similar to our 
local strategy. This is consistent with experimental work [33] (see also [19]) showing 
a dip in the fixed-nitrogen level approximately halfway between two widely separated 
heterocysts. 

Qualitatively our best strategy is easy to state: "quickly differentiate cells that are 
locally starving of fixed-nitrogen" . Effectively implementing the strategy is not trivial. 
We note two complications. The first is that our local strategy works for ongoing 
differentiation of heterocysts, during steady-state growth. In that regime, long after 
efN deprivation, a dedicated mechanism to avoid Mch is not needed — because nearby 
cells begin to starve at different times. However, in the first burst of differentiation we 
need to avoid Mch with a temporary mechanism to prevent differentiation of the entire 
filament. Indeed, the initial heterocyst spacing is distinct from the steady-state pattern 
[8] and cyanobacterial filaments are thought to use diffusible inhibitors derived from 
Pats to suppress Mch [7|. The second is that we have simplified heterocyst commitment 
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[S] and subsequent delay before du novo nitrogen fixation into a single delay r. While 
we find that a smaller r always leads to faster filament growth, heterocysts do take 
significant time to both commit and to begin to fix nitrogen [3l El ETJ [38] (though 
see |36j). While we have explored a correspondingly large range of r G {1 — 20}h, 
with similar results throughout, the details may be expected to change with a more 
detailed model of heterocyst commitment. However our understanding of heterocyst 
commitment timing with respect to local fN concentrations remains crude, and we feel 
that a unified delay r is an appropriate simplification. 

We predict a plastic developmental response of heterocyst frequency to levels of 
efN, where the heterocyst frequency / decreases rapidly with pefN- This implicitly 
assumes that there is no threshold extracellular fixed-nitrogen concentration above 
which heterocysts will not differentiate, and more broadly that there is not a fixed 
developmental pattern of heterocysts. Supporting a plastic developmental response, 
early work by Fogg [13] showed a time- dependent heterocyst frequency that increased 
with decreasing levels of efN, Ogawa and Carr [30] found heterocyst frequency and 
nitrate concentration to be inversely related, and Horne et al [31] found the same 
correlation in field studies. Fogg also observed significant heterocyst development at 
40/iM extracellular fixed-nitrogen (Fig. 2 of [13]), though not in steady-state conditions. 
Plastic developmental responses are observed in systems ranging from plants [40j to 
the brain [H]. We feel that filamentous cyanobacteria are particularly amenable to 
exploring and understanding the adaptiveness of this plastic response. The time is now 
ripe to undertake constant external fixed-nitrogen concentrations with single-filament 
studies in the chemostat-like environment of microfiuidic systems, such as filamentous 
cyanobacteria within the long channels used in studies of persister cells of E. coli [421143] . 

We note that the heterocyst spacing patterns themselves are not adaptive. Rather 
the patterns reflect an adaptive local strategy for placing heterocysts close to locally 
starving cells. This starvation is a combination of distance from heterocysts and 
fast growing cells. Local placement is particularly adaptive when fN leaks from the 
cyanobacterial filament (see [13], [111 [151 [161 [TTl HH] , ^ discussed in the introduction). 

Our heterocyst placement strategy is deliberately simple, and does not include any 
diffusible inhibitors such as peptides derived from PatS or HetN [Tll44l[45]. We anticipate 
that the small secondary peak that we observe at smaller separations in Figs. [6] (a)-(d) 
may be suppressed in vivo by the action of such diffusible inhibitors. Since HetN appears 
to act after initial heterocyst development [10] , it may indeed be why the secondary peak 
is not actually observed [8]. However, we note that the separation corresponding to the 
secondary peak increases with the average separation, so that the secondary peak may 
emerge in experimental studies at larger levels of Pe/Af if inhibition due to PatS and 
HetN has a fixed range. 

Our model and its results assume that, apart from fixed-nitrogen, there are no 
nutrient or metabolite requirements that limit growth differentially along the filament. 
Doubtless, there are experimental conditions where this assumption is not warranted. 
Certainly, there are interesting nutrients that we have not included in our model. For 
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example, the supply of membrane potential and/or of carbohydrates from vegetative 
cells to heterocysts ^61 W7\ was not considered. Nevertheless, while carbohydrate 
transport would be expected to starve the interior cells of clusters of contiguous 
heterocysts and reduce their capacity for nitrogen-fixation, our primary focus was on 
the local heterocyst placement strategy under conditions of non-zero external fixed- 
nitrogen — where no clusters of heterocysts are observed and heterocysts are broadly 
spaced. However, carbohydrate limitation would affect random placement strategies, 
which might then have even more reduced growth compared to our best, local strategy. 

For cyanobacterial filaments implementing a local heterocyst placement strategy we 
have several experimentally testable predictions that we expect to be observed in steady- 
state conditions — independent of the specific parameterization of our quantitative 
model. To explore them in the lab, external concentrations of fixed- nitrogen (pe/Af) 
would need to be controlled by microfluidic devices [l2l |13] or flow cells. The first is 
that heterocyst frequency /* will decrease rapidly and continuously as pe/w is increased 
from zero. As shown in Figs. H] (b) and (d), the dependence is markedly non-linear with 
larger leakage rates of fN from the growing filament. The second, related, prediction is 
that the peak of the heterocyst spacing distribution will increase as pefN decreases. This 
is illustrated in Figs. [6] (a)-(c). We also expect that the width, or standard deviation, 
of the heterocyst spacing will increase as the mean spacing increases. Indeed, for our 
model, the standard deviation and mean spacing are approximately proportional. 

Our results may also have implications for quantitative models of biogeochemical 
cycling of fixed-nitrogen in marine and lake environments, where filamentous 
cyanobacteria can play a significant role. We note that a number of existing models 
[29l l48l |49] use growth rates and nitrogen-fixation rates that are independent of the 
biologically available fixed-nitrogen. We find that, with a local heterocyst placement 
strategy, the filament growth rate fi* is approximately (but not exactly) independent 
of PefN — see Fig. H] (a) and (c), consistent with the approximately constant growth 
reported for chemostat studies by Elder and Parker [32]. However we expect that du 
nuovo nitrogen- fixation to be proportional to /*, which in turn strongly decreases with 
increasing Pe/N- 
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